Clustering in disordered ferromagnets: The Curie temperature in diluted magnetic 

semiconductors 
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We theoretically investigate impurity correlation and magnetic clustering effects on the long-range 
ferromagnetic ordering in diluted magnetic semiconductors, such as Gai-^Mnz As, using analytical 
arguments and direct Monte Carlo simulations. We obtain an analytic formula for the ferromagnetic 
transition temperature T c which becomes asymptotically exact in the strongly disordered, highly 
dilute (i.e. small a;) regime. We establish that impurity correlations have only small effects on T c 
with the neutrally correlated random disorder producing the nominally highest T c . We find that the 
ferromagnetic order is approached from the high temperature paramagnetic side through a random 
magnetic clustering phenomenon consistent with the percolation transition scenario. 
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MOTIVATION AND INTRODUCTION 

The effect of disorder on magnetism is an old and im- 
portant problem in condensed matter physics and ma- 
terials science. Worldwide current activity 1-3 in diluted 
magnetic semiconductors (DMS) where cationic substi- 
tutional doping (by a few percent) of a semiconduc- 
tor with magnetic impurities (e.g. Gai-^Mn^As with 
x ~ 0.01 — 0.1) seemingly leads to an intrinsic ferromag- 
netic material with a relatively high Curie temperature 
(T c w 1QQK — 200-fiQ, has renewed vigorous recent inter- 
est in this important subject. In particular, the intrinsic 
mechanism of ferromagnetism (i.e. the physics underly- 
ing the long-range ordering of the randomly distributed 
magnetic impurity moments) as well as the precise role 
of disorder (e.g. arising from the random distribution of 
magnetic Mn ions at the cation sites in the GaAs lat- 
tice) are substantial questions, both from the fundamen- 
tal perspective of understanding the competition between 
disorder and magnetic interactions as well as the techno- 
logical perspective of the fledgling subject of 'spintronics' 
(or spin electronics) where the projected seamless inte- 
gration of magnetics and electronics on a single ferromag- 
netic semiconductor chip could lead to a new paradigm in 
microelectronics 4 . In this article, we study theoretically 
the competition between quenched disorder (i.e. the ran- 
dom spatial distribution of the magnetic impurities) and 
magnetic correlations in DMS materials using analyti- 
cal arguments supported by extensive large-scale Monte 
Carlo simulations on a disordered Hciscnberg model. Our 
main findings are that (1) site disorder by itself has a rel- 
atively small effect on the magnitude of the DMS ferro- 
magnetic transition temperature T c , and (2) global ferro- 
magnetic ordering in DMS is approached from the high- 
temperature (i.e. T > T c ) paramagnetic phase through 
the formation of random disconnected magnetic clusters 
of increasing size as T approaches T c from above, which 
coalesce at T = T c leading to a magnetic percolation 
transition, similar to what was postulated in ref. 5 . We 
believe that our theoretical findings are quite general, and 



should apply to all diluted site disordered ferromagnetic 
materials including doped magnetic oxides 6 ' 7 , which are 
insulators or poor metals (with short mean free paths 
comparable to lattice spacings). In particular, we spec- 
ulate that the results presented herein are sufficiently 
general to be valid for all site disordered ferromagnets 
(with localized impurities providing the underlying mag- 
netic moments) in the dilute (or, equivalently, strongly 
disordered) regime, though we discuss our work entirely 
in the context of Gai^Mn^As since this is by far the 
most extensively studied DMS material experimentally 
and theoretically. 

A key feature 1-7 of DMS materials is that the fer- 
romagnetic interaction between the impurity magnetic 
moments is an indirect exchange interaction medi- 
ated by the semiconductor carriers (which are holes in 
Gai-zMn^As). The 'standard' model for DMS ferromag- 
netism was developed 8 more than forty years ago, and has 
recently been rediscovered 9 . This model, which we refer 
to as the VCA model, is a simple mean-field model in 
both the spatial disorder (i.e. a continuum virtual crystal 
approximation, VCA) and the thermal fluctuations (i.e. 
a Curie- Weiss mean field theory for the ferromagnetic 
ordering of the impurity moments) with the magnetic 
interaction between the impurity moments being medi- 
ated by the indirect exchange interaction arising from the 
carrier-induced RKKY-Zener mechanism 1-3 . Although 
the extreme simplicity of the VCA model, for example 
its analytic prediction of T c cx mJ^ d n\^x (where J p d is 
the effective p-d exchange interaction strength between 
the p-type holes in GaAs and the Mn d-level, and m, n c 
are the effective hole mass and density, respectively) , has 
led to its widespread application in the DMS literature, 
we emphasize here that the model is conceptually wrong 
since it predicts a monotonically increasing ferromagnetic 
transition temperature with increasing carrier density, 
impurity density, and magnetic coupling strength - in 
fact, all three trends are conceptually incorrect 10 (and 
the simple VCA model provides no mechanism for esti- 
mating its regime of validity in the J p( i, n c , x, T, and 
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disorder parameter space). The questions (e.g. impurity 
correlation, magnetic clustering) being asked in this work 
are simply beyond the scope of VCA due to its complete 
neglect of site disorder and thermal fluctuations. We also 
find that in the strong disorder limit (i.e. highly diluted 
case of x <C 1), where the VCA fails most miserably, 
T c (J p d, x,n c ) has a qualitatively different behavior than 
that predicted by the standard DMS VCA mean field 
theory. 

ASPECTS OF THE CALCULATION 

We assume Mn ions randomly occupy only Ga sites in 
the GaAs zinc-blende lattice (fee) with a lattice constant 
a. We operate in the (weakly) metallic limit and assume 
the carrier-mediated effective Mn-Mn indirect exchange 
interaction to be of the RKKY form. However, with DMS 
systems being at best poor or bad metals with a mean 
free path typically of the order of l-2Awhich is less than 
a, it is important to include the effects of a finite carrier 
mean free path; we do this by introducing a cutoff Z(< a) 
in the range of the RKKY function. The large value (S — 
5/2) of the impurity moment spins allows a treatment of 
the Mn spins as classical Heisenberg spins. Hence, our 
Hamiltonian is given by 

H = ^J(r ij )S i -S j , (1) 

where is the separation between moments i and j, and 
J(r) is the damped RKKY range function given in three 
dimensions by 

J(r) = J Q e- r/l r- 4 [sm(2k F r) - 2k F r cos(2fc F r)]; (2) 

A;f = (|7r 2 n c ) 1 / 3 is the Fermi wave number and n c is the 
hole density. Jo(> 0) is related to the local Zcner cou- 
pling J p d between the Mn local moments and hole spins, 
J cx m Jp d where m is the hole effective mass (for conve- 
nience, spin units are used with the factor S 2 = (5/2) 2 
absorbed into Jo). While I sets the range of interac- 
tion between moments, l s = n i (where rii is the vol- 
ume Mn density) gives the mean inter-impurity separa- 
tion and provides a measure of the disorder strength with 
l s w a for weak disorder and l s » a for strong disorder. 
Note that by definition the strong disorder limit is equiva- 
lent to the dilute limit, which is the situation appropriate 
for DMS materials. For weak disorder with a long range 
interaction (i.e for I 3> {l s ,a}), Curie Wiess Mean Field 
Theory (MFT) fares well (though lattice MFT 11 is an 
important refinement in addressing the discrete crystal 
structure); we show here that strongly disordered ferro- 
magnets can also be viewed in terms of a limiting the- 
ory in which the Curie temperature asymptotically ap- 
proaches 

T c = ^Lj(2r c n; 1/d ), (3) 



where rj is a constant of order unity and r c is the di- 
mcnsionless critical radius for random sphere percolation 
with r c = .4436 12 (We concentrate entirely on the cur- 
rently experimentally relevant three-dimensional (d = 3) 
DMS systems in this paper). It is clear that this T c ex- 
pression, T c ~ J^n// 1 ^), bears little resemblance to the 
VCA result T c cx n, J J(r)d d r. 

To calculate quantities such as T c and the spontaneous 
magnetization M, we use a hybrid Monte Carlo technique 
combining Wolff Cluster 13 and Heat Bath 14 moves. The 
Wolff moves overcome critical slowing down near T c , and 
Heat Bath moves ensure that ergodicity is attained ef- 
ficiently by thermalizing spins weakly coupled to their 
neighbors. We calculate T c by locating the crossing tem- 
perature of Binder cumulants C/ 15 for two different sys- 
tem sizes; in our case, U = 1- (1/3)[(|M| 4 )]/[(|M| 2 )] 2 , 
where angular and square brackets denote thermal and 
disorder averaging, respectively. Since we invariably have 
(JVi mp ) = riiL 3 > 2 x 10 4 , we have found it useful to ex- 
ploit our large system sizes; calculating magnetization 
moments |M| l by simply raising the total magnetiza- 
tion M to the ith power disposes of information that 
can be gleaned by dividing the system into n s subsys- 
tems and averaging |M| J evaluated for each of these. 
(n s = (iVimp} 1 / 2 has proven to be a reasonable choice). 
Operating in this manner, we use as few as 50 total 
Monte Carlo sweeps (n cq — 25 for equilibration and 
n av = 25 for statistical sampling) with no compromise 
in accuracy. We have checked for each calculation the 
adequacy of the equilibration stage by comparing "cold 
starts" (systems prepared with perfect ferromagnetic or- 
der) and "hot starts" (systems with initially randomized 
spins), finding in each case agreement to within random 
statistical error (confined to less than 1% for all results 
reported here). 

CORRELATED MAGNETIC MOMENTS 

Relaxing the requirement that magnetic impurities re- 
side on the discrete fee lattice in favor of continuously 
distributed Mn is a simple way to examine the sensitiv- 
ity of the ferromagnetic state to how the Mn ions are dis- 
tributed. In principle, the discrete and continuum cases 
represent very different situations. Even so, it its evi- 
dent in Figure 1 that Curie temperatures for the discrete 
(T lat ) and continuum (T^° nt ) cases (shown for more than 
a decade of x values) are in close agreement. The main 
graph shows critical temperatures directly for l/a = 0.5, 
inset (a) depicts the ratio of continuum and lattice T c 's 
for l/a — 0.5, and inset (b) displays the same result for 
the even smaller cutoff l/a = 0.25. In the continuum 
model, it is possible for magnetic impurities to cluster 
more tightly than in the discrete case, and the spins in 
such clusters are therefore more strongly aligned. How- 
ever, the enhanced intra-cluster coupling of the compact 
clusters tends to be offset by a weaker coupling to neigh- 
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boring spins (since the overall constancy of x requires 
that compact structures be more isolated, on average, 
than more diffuse Mn assemblies). The close agreement 
of T c cont and T 1 ^ in Figure 1 suggests that these two ef- 
fects nearly precisely cancel. It is important to note that 
the close similarity of continuum and discrete T c 's does 
not justify the VCA result 9 . The latter, given by 



T vcA = ijT Jo x 



! J(r)dr 



(4) 



does not take into account thermal and disorder fluc- 
tuations (important in the dilute DMS regime). The 
disparity between the VCA and Monte Carlo results is 
highlighted in panel (c) of Figure 1, and is even greater 
if the damping factor e~ T l l is not taken into account. 
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FIG. 1. In the main graph, T c cont (continuum) and 
Tl at (discrete) are graphed with continuum T c 's shown as 
closed circles and discrete T c 's plotted as open circles; 
nc/rii = 0.1 and l/a = 0.5. In inset (a), ratios T c cont /T] at 
are shown; again, n c /rii =0.1 and l/a — 0.5. The graph in 
inset (b) shows T c cont /T c lat for n c /n t = 0.1 and l/a = 0.25 with 
the horizontal dark reference line indicating unity. Panel (c) 
shows VCA (solid and broken lines for damped and undamped 
cases, respectively) and Monte Carlo T c 's on the same graph 
with n c /rii — 0.1 and l/a = 1.0 (the damping refers to the 
inclusion of mean free path effects). 



impurity clustering 

To examine impurity correlations in a way that is con- 
tinuously tunable, we introduce a local attractive in- 
teraction between impurity moments yielding the con- 
figurational energy -E con f = — Eq X^7in(b° n ds) wnere 
^2 nn (bonds) is the total number of bonds shared by 
neighboring Mn ions. We prevent the formation of a 
single large cluster by limiting the lateral dimension of 
clusters, requiring in this case that the size along any 
axis be smaller than a. We vary the extent of impu- 
rity clustering via an annealing temperature Ta; large 
Ta values lead to weaker clustering, while in the low Ta 
limit, impurities tend to be more exclusively bound in the 
complexes mentioned above. Configurations are prepared 
by subjecting initially neutrally correlated samples to an 
equilibration stage at Ta (consisting of 500 Metropolis 
sweeps), and the extent of clustering in the resulting 
sample is characterized by (n s i ze ), the mean cluster size. 
Impurity moments locked in clusters tend to align due 
to their proximity, but the large inter-cluster spacings 
makes the alignment of cluster moments with each other 
more readily disrupted by thermal fluctuations. In fact 
as x is decreased, one ultimately obtains an essentially 
paramagnetic state composed of large but very weakly 
coupled cluster moments. In Fig. 2(a) and Fig. 2(b) 
the impact of clustering on T c is shown for two differ- 
ent Mn concentrations x for clustering ranging from very 
weak to quite strong. For x = 0.03 (the case depicted in 
Fig. 2(a)), the overall Mn density is sufficiently low that 
one does see a slight dimunition of T c for increasing clus- 
tering. However, when x is doubled to x = 0.06 (shown in 
Fig. 2(b)), the T c curve is essentially flat. Hence, even for 
fairly small Mn densities, the effects of enhanced intra- 
cluster alignment and diminished inter-cluster coupling 
are roughly the same in magnitude, leading to a weak net 
effect. Thus impurity correlation or clustering does not 
seem to have much effect on T c . We mention as an impor- 
tant caveat that in Gai-^Mn^As, two closely placed Mn 
ions (e.g. a substitutional and interstitial Mn atom in the 
same unit cell) experience very strong nearest-neighbor 
short-range direct antiferromagnetic exchange coupling. 
This would lead to 11 a strong suppression of ferromag- 
netic T c at large values of x if there is substantial Mn 
clustering. This effect is not included in our model. 
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FIG. 2. Curie Temperatures for impurity clusters (panels 

(a) and (b)) and impurity superlattices (panel (c)). In (a) and 

(b) , where x is respectively 3% and 6%, T c is plotted versus 
(^size), the mean cluster size. In (c), T c is graphed versus 
(rd), the mean deviation from superlattice sites. For (a), (b), 
and (c), n c /n; = 0.1 and l/a — 0.5. 

impurity superlattices 

In a sense opposite to inducing impurity clustering 
is arranging Mn ions in superlattices, which tends to 
break up impurity clumps to form a homogeneous lattice. 
Again, we begin with perfect superlattices and introduce 
deviations by annealing at the "disordering temperature" 
Td; s ; the mean displacement of impurities (rd) from their 
sites on the superlattice characterizes the resulting con- 
figuration (very small values of (rd) correspond to nearly 
perfect superlattices while (rd) ~ a indicates an essen- 
tially "melted" state with randomized impurities). Curie 
temperatures for a broad range of (rd) are shown in panel 

(c) of Fig. 2 (our superlattice is simple cubic with lattice 
constant 2a, which corresponds to x = 0.031). While the 
calculated T c in (c) slightly increases as the superlattice is 
disrupted, the effect is relatively modest, consistent with 
previous work comparing random disorder and pure su- 
perlattice configurations 17,18 . The results suggest that 
while an excessive degree of clustering can lower T c , the 
absence of any inhomogeneities in the impurity configu- 
rations is also deleterious to the ferromagnetic state. As 
the transition is made from a perfect superlattice to more 
disordered configurations, more and more impurities are 
displaced and can act as a bridge between neighboring su- 
perlattice sites, thereby increasing the effective coupling 
between neighboring Mn ions still occupying superlattice 
sites. Ultimately, the T c curve saturates for large (rd) 
where superlattice modulations vanish. It thus appears 
that the uncorrelated impurity distributions occupy a 
"happy" medium with a critical temperature higher than 
that of either strongly clustered configurations or states 
with superlattice structure. Thus random disorder mod- 
estly enhances DMS T c . 

THE STRONG DISORDER LIMIT 

Finally, we examine the regime of very strong disor- 
der (i.e. l s 3> {I, a}), finding that the Curie temperature 



exhibits well-defined asymptotic behavior in the large l s 
limit corresponding to the percolation of randomly placed 
spheres. We first exploit the strong similarity of the be- 
havior of discrete and continuum systems shown in Fig- 
ure 1 by relaxing the discrete occupancy requirement, 
although this condition is still imposed in the numerical 
calculations. For convenience, we rescale the linear di- 
mensions in such a way that l s is equal to unity, thereby 
mapping our problem to a new ferromagnetic system 
where the range of J(r ) is shortened to I = l/l s . Hence, 
for large l s , J(r ) varies quite rapidly and as a result the 
coupling between pairs of spins tends either to be very 
large or very small. The strong disorder limit thus can be 
viewed as the regime in which spins even slightly closer 
than the "correlation radius" r corr are strongly correlated 
whereas for separations greater than r corr , the coupling 
rapidly becomes small relative to fc^T. This is essentially 
a real space renormalization group (RG) argument for the 
scaling of the spin coupling. This observation motivates a 
comparison with the problem of randomly placed spheres 
(distributed with unit density n — 1), which form perco- 
lating networks when their radii exceed r c ; it is natural 
to identify the critical value of r corr with 2r c and to write 
fcsT* = J(r corr ) for the ferromagnetic transition temper- 
ature 19 . However, a constant i] of order unity is needed to 
take thermal fluctuations into account 20 and one actually 
has T c = -P-J(2r c l s ), where the final form is expressed in 
terms of the original, unmodified units. Our MC simula- 
tions give essentially exact (within 1 — 2%) quantitative 
agreement (shown in Fig. 3 for the x = — 0.1 range 
we checked) between the numerical MC T c results and 
the theoretical T* with both showing essentially singular 
x 4 ^ 3 exp(— en" 1 ' 3 ) behavior for small x. 




FIG. 3. The graph shows together the numerical MC T c 
results (open circles) and the the theoretical T* (solid line) 
for n c /rii = 0.1 and l/a — 0.5 

We show explicitly in Fig. 4 that the approach to T c 
from above is marked by the formation of magnetic clus- 
ters of correlated spins by visualizing these clusters us- 
ing the Swendsen-Wang algorithm 21 (which partitions 
the entire system into spin clusters). To facilitate view- 
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ing, we examine a 2D system, a dilute site-disordered 
Ising model, where (x = 0.05 and l/a = 1.0). In Fig. 4, 
Swendson-Wang snapshots are depicted for four temper- 
atures. One can see the merging of disconnected clus- 
ters to form a spanning cluster at T c , though substantial 
clustering occurs even well above T c . This is very simi- 
lar to the magnetic polaron percolation picture of DMS 
ferromagnetism theoretically envisioned in Ref. 5 and the 
qualitative discussion in Refs. 17 ' 18 . We point out, how- 
ever, that the T c formula in our RG argument (validated 
by our MC simulations) differs somewhat from the corre- 
sponding percolation theory result 5 . We emphasize that 
although our discussion has been in the DMS context, 
our arguments are not dependent on particular traits of 
the RKKY range function. 




(c) T = T c (d) T = .5T C 

FIG. 4. The four panels (a,b,c, and d) display clusters 
of correlated spins for x = 0.05, J(r) = Joe~ r /' where 
l/a — 1.0. Each panel corresponds to a different tempera- 
ture with T/T c = {4.0,2.0,1.0,0.5} for (a), (b), (c), and (d) 
respectively. Dark/open circles represent "up" /"down" spins. 

CONCLUSIONS 

In conclusion, we have examined the effect of impu- 
rity correlation in diluted magnetic semiconductors. We 
hnd that the presence of correlations among impurity 
positions affects the critical temperature T c only mildly, 
with the calculated T c being optimal for neutral impurity 
correlations (i.e. random disorder). In addition, we have 
verified via explicit numerical calculation that just as fer- 
romagnetic behavior in the weak disorder mean field limit 
is amenable to mean field approaches, the strong disor- 
der regime also can be understood in terms of a limiting 



theory, random sphere percolation. We develop a new 
theory for the dependence of T c on the magnetic impu- 
rity density by mapping the discrete random problem to a 
continuum problem using real space RG arguments, and 
validate the theory by Monte Carlo simulations. This 
work is supported by US-ONR. 
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